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Abstract 

By employing the N-barrier method developed in the paper, we establish a new 
N-barrier maximum principle for diffusive Lotka-Volterra systems of two compet¬ 
ing species. As an application of this maximum principle, we show under certain 
conditions, the existence and nonexistence of traveling waves solutions for systems 
of three competing species. In addition, new (1,0, 0)-(m*, u*, 0) waves are given in 
terms of the tanh function provided that the parameters satisfy certain conditions. 


1 Introduction 


Species diversity refers to the number of different species and abundance of each species 
that live within an ecological system. To be more specific, species diversity takes into 
consideration species richness and species evenness; the former is defined as the total 
number of different species and the latter the variation of abundance in individuals per 
species. The importance of species diversity to an ecological system lies in the fact that 
the ecological system with a greater species diversity may function more efficiently and 
productively since more resources available for other species within the ecological system 
will be made. Therefore, the study of species diversity has been extensively investigated 
via both field research and theoretical approaches. 

As a suggestive example of species diversity, we investigate the situation where one 
exotic competing species (say, W) invades the ecological system of two native species 
(say, U and V) that are competing in the absence of W. Then a problem related to 
competitive exclusion (0, mi, H, izDi. izi, H) or competitor-mediated coexistence 
(0, 1231, H) arises, and a mathematical model for this situation can be provided by 
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the diffusive Lotka-Volterra system of three competing species (n, im, n. iia, la, 

m, la. ra. ffl. ra. pd 

' Ut = di Uyy + M (cTi — Cii M — Ci2 U — Ci3 ta), 2 / € t > 0 , 

< Vt = d2Vyy + V {a2 - C21U - C22V - C22,w), e M, t > 0, (1.1) 

= da + w (0-3 - C31M - C32U - C33 w), y t > 0 , 


where u{y,t), v{y,t) and w{y,t) denote the population densities of U, V and W at time 
t and position y. The parameters di, ai, ca {i = 1, 2, 3), and Cij {i,j = 1, 2, 3 with i 7 ^ j), 
which are all positive constants, stand for the diffusion rates, intrinsic growth rates, intra¬ 
specific competition rates, and inter-specific competition rates, respectively. In particular, 
the result in [7j indicates from the viewpoint of competitor-mediated coexistence, the 
coexistence of strongly competing species in the presence of an exotic competing species. 


We begin with a two-species system of (1.1) in the absence of W, i.e. 


Ut = di Uyy + u{ai — cii M — C 12 n), 2 / G K, f > 0 , 

Vt = d 2 Vyy + v{a 2 - C21 U- C22V), 2/ e M, f > 0 . 

Imposing the zero Neumann boundary conditions 

du dv 

— = — = 0 , X edVt 

ov ov 


(1.2) 


(1.3) 


and suitable initial conditions 


u{x,Q) = u{x), n(a:, 0 ) = v{x), x E fl 


(1.4) 


on ( 1 . 2 ) with the entire space M is replaced by a bounded and convex domain fl, we 
conclude from na and m that any positive solution {u{x,t),v{x,t)) of such a initial¬ 
boundary value problem converges to either (^,0) or (0, ^) when U and V are strongly 
competing, i.e. when the following condition hold: 


(Ti (72 (72 (7i 

Cll C21 C22 C12 


(1.5) 


In this case. Cause’s principle of competitive exclusion occurs between the two species U 
and V when [/ and V competing for the same limited resources cannot stably coexist; 
one will prevail and the other is excluded. When the influence of diffusion in (1.2) is 
disregarded, ( 1 . 2 ) becomes 


{ Ut = u (ai - cnu - Ci 2 v), t > 0, 

Vt = V ((72 - C21 M - C22 v), t > 0. 


( 1 . 6 ) 


It is readily seen that (1.6) has four equilibria: ei = (0,0), 62 = (^,0), 63 = (O, 


and e 4 = {u*,v*), where {u*,v*) = ( 


C22 (Tl—C 12 «T2 Cll «T2—C21 tJi 
Cll C22—C12 C21 ’ Cll C22—C12 C21 


) is the intersection of the 
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two straight lines cxi — Cn m — C12 n = 0 and ct2 — C21 n — C22 n = 0 , whenever it exists. In 
the diffusion-free case, we can classify the asymptotic behaviour of solutions {u{t),v{t)) 
for (1.6) as f —)■ cxD depending on and Cij {i,j = 1,2), as described in: 


Proposition 1.1 ([S]). Suppose that {u{t),v{t)) is a solution of (1.6) with initial condi¬ 
tions M(0),n(0) > 0. We have 


(i) (Competitive exclusion) ^ ^ ^ then \im{u{t),v{t)) = (^,0); 

(a) (Competitive exclusion) ^ ^ o,nd ^ then \im{u{t),v{t)) = (O, 


(in) (Strong competition) ^ ^ o,nd ^ then \im{u{t),v{t)) = (^,0) or 

( 0 , depending on the initial condition; 


(iv) (Weak competition) ^ ^ ^ then \im{u{t),v{t)) = {u*,v*). 


As in case {Hi) mentioned previously. Cause’s principle of competitive exclusion also 
occurs in cases {i) and {ii): one species always wins and the other species become extinct 
in the long run. It is easy to see that we do not need to treat one of cases (z) and {ii) 
since one of the two cases is obtained from the other by exchanging the roles of u and v. 
For the case of weak competition, the Lotka-Volterra model (1.6) predicts that the stable 
coexistence state {u*,v*) exists only when intra-specific competition has a greater effect 
than inter-specific competition. 

We shall assume throughout, unless otherwise stated, that either strong or weak 
competition occurs between the two species U and V: 


• (Strong competition) [SI: ^ ^ and — > —; 

• (Weak competition) [Wl: ^ ^ and ^ 

^ Z L J dl C 21 C22 C 12 

in investigating traveling wave solutions 


= {u{x),v{x)), x = y-0t 


(1.7) 


of (1.2). Here 9 is the propagation speed of the traveling wave. We note that u*,v* > 0 
if and only if either [S] or [W] holds. Substituting (1.7) into (1.2) yields the following 
system of ordinary differential equations 


{ di Uxx + 9ux + u (cTi — cii M — C 12 n) = 0, a: G M, 

d2 Vxx +9Vx-kV ((72 - C 21 M - C 22 t") = 0, X G M. 


( 1 . 8 ) 


The problem as to which species will survive in a competitive system is of importance 
in ecology. In order to tackle this problem, we use traveling wave solutions of the form 

(PI). 
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In this paper, we treat the following bonndary valne problem 


" dl Uxx + 9 Ux + u (cTi — Cii M — Ci2 v) = 0, x G M, 

< d2Vxx + Ovx + V {a2 - C21U - C22V) = 0 , X e M, (1.9) 

^ (m, n)(-cx)) = 62, (m, n)(+oo) = 64. 


and call a solntion {u{x),v{x)) of (1.9) an ( 62 , e 4 )-wave. 


Under various parameters, monotone ( 62 , e 4 )-waves are investigated via different ap¬ 
proaches (for instance, [11], |2I], [22]). It is not clear from the assumptions of parameters 
given in the references above that the monotone ( 62 , e 4 )-waves have the property u* > v* 
or not. Let us see what happens when the answer is affirmative. If u* > v* holds, then we 
easily see that, since v is monotonically increasing and u is monotonically decreasing the 
profile of u lies completely above the profile of v. Accordingly, in this case u dominates the 
entire habitat M. Indeed, it has been proved that there exist two types of ( 02 , e 4 )-waves, 
one with u* > v* and the other with u* < v*, by giving exact ( 02 ,e 4 )-waves in [H]. We 
note that in the case u* > v*, the phenomenon exhibited by the dominance of u in the 
entire habitat M is unique and is yet to be explored. In particular, exact (ei, e 2 )-waves, 
( 01 , 64 )-waves, ( 62 , 63 )-waves and ( 02 , 64 )-waves are given in [19] under certain conditions 
on the parameters by applying judicious ansatzes for exact solutions. 

When the inhabitant of the two competing species U and V is resource-limited, the 
investigation of the total mass or the total density of the two species U and V is essential. 


This gives rises to the problem as to the estimate of the total density u(x) + v(x) in (1.9). 
In addition, another issue which motivates us to study the estimate of u(x) + v(x) is the 


measurement of the species evenness index J for (1.9). As J is defined via the Shannon’s 
diversity index "H ([2], [H], [2T], [22]), i.e. 


J = 


n 

ln(s)’ 


( 1 . 10 ) 


where 


91 = - ^Pi ■ In(pi), 


( 1 . 11 ) 


i=l 


s is the total number of species, and pi is the proportion of the Ath species determined 
by dividing the number of the Ath species species by the total number of all species, the 


species evenness index J' for (1.9) is given by 


J = 


(ln2)('U -f v) 


u In 


u 


U + V 


n In 


V 


u + V 


( 1 . 12 ) 


One of our primary goals in this paper is to address the problem of giving a priori 


estimates of u{x) + v{x), which is involved in the calculation of JL in (1.12). On the 


other hand, we also are concerned with the following question when a priori estimates of 
u{x) + v{x) are given: 


Q1 


: How does the estimate of u{x) -f r;(x) depend on the parameters in (1.9)^ 
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In [U], upper and lower bounds of u{x) + v{x) are given when the two diffusion rates 
di and d 2 are equal. However, the approach employed in [ 6 ] to obtain upper or lower 
bounds for M + n cannot be applied directly to the case where the diffusion rates di and 
d 2 are not equal. 


Q2; In when di ^ d 2 , 


can upper and lower bounds of u + v be obtained? 


As for the answer to Q2, it seems as far as we know, not available in the literature. To 
give an affirmative answer to this question, we develop a new but elementary approach. 
Moreover, employing this approach leads to a affirmative answer to the following question 
which is more general than Q2: 


Q3: In ( |1.9[ ), when di ^ d 2 , 
are arbitrary constants, be given 


can upper and lower bounds of au + (3v, where a, ft > 0 
Q 


By adding the two equations in (1.9), we obtain an equation involving p{x) = a u{x) + 


ft v{x) and q{x) = dia u{x) + d 2 ftv{x), i.e. 

0 = a [diUj;,, + 6 Uj; + u (di -ciiu- v)) + ft (^2 + Ov,, + v (a 2 - C 2 iU- C 22 v)) 

= q"{x) + 6p'{x) + au{ai - cnu - Cu v) + ftv{a 2 - C 21 u - C 22 v) 


= q"{x) + 6p'{x) + F{u, v). 


(1.13) 


For the sake of clear exposition we shall assume from now on that F{u,v) = au {ai — 
Ciiu — Ci 2 v) + ft V ((72 — C 21 M — C 22 V). The case where di = ^2 or p{x) is a non-zero 
constant multiple of q{x) has been considered in [B]. A mathematical difficulty arises as 
a consequence of the fact that the approach used in [ 6 ] cannot be directly applied, when 
di 7 ^ d 2 in the last equation since p{x) no long can be written as a constant multiple of 
q{x) and such p[x) and q{x) are involved in a single equation like (1.13). The approach 


proposed here can be employed to give estimates of pfx) in the case where p{x) and q{x) 


are involved in the single equation (1.13). One of the main results of this paper is the 


N-barrier maximum principle (Theorem 2.1), which gives an affirmative answer to Q3. 
The rest of the paper is organized as follows. In the next section, the main results of 


this paper, including the N-barrier maximum principle (Theorem 2.1) and two applica¬ 
tions (Theorem 2.2 and Theorem |2.3[ ) to the system of three species (1.1), are summarized. 
We prove in Section the N-barrier maximum principle by means of the construction of 
N-barriers depending on various conditions. Under certain restrictions on the parameters, 
the existence of exact traveling waves solutions for (1.1) (Theorem |4.1 ) is presented in 
Section]^ Finally, Sectionis devoted to the proofs of Theorem 2^ and Theorem 2^ 


2 Statement of main theorems 


Theorem 2.1 (N-barrier maximum principle). Under either [S] or [W], we assume 
that {u{x),v{x)) with u{x),v{x) > 0 forx G M satisfies the boundary value problem (1.9). 
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For any a, (3 > 0, we have 


< au{x) + (3v{x) < q*, a; G M, 


where 


and 




■ (7i (72 ■ 

, (3 min 

■ (72 (7i ■ 


' di (^2 

g* = min 

a min 

■Cll ’ C 21 . 

-C 22 ’ C 12 - 

min 

.d2 di- 




■ (7i (72 ■ 

, (3 max 

■ (72 (7i ■ 


d\ 6^2 

max 

a max 

■Cll ’ C 21 . 

-C 22 ’ C 12 - 

max 

- G?2 d\ - 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


In particular, 


(I) (Lower bounds for q = q{x)) when the differential equations in (1.9) are replaced 
by the differential inequalities 


di Uxx + 9 Ux + u {ai — Ciiu — Ci2 1 ") < 0 , x G M, 
d2 Vxx +9vx + v ((72 - C 21 U- C22 v) < 0 , X G M, 
we have au{x) + (3v{x) > g* for x G M; 


(2.4) 


(II) (Upper bounds for q = g(x )) when the differential equations in (1.9) are replaced 
by the differential inequalities 


di Uxx + 9 Ux + u {ai — Ciiu — C 12 x) > 0, x G M, 

d 2 Vxx +9vx-\-v ((72 - C21U- C22 v) > 0 , X G M, 

we have a u{x) + (3 v{x) < q* for x G M. 


(2.5) 


We would like to add a few remarks concerning Theorem 2.1 

(i) In addition to the diffusion rates di and d 2 , upper and lower bounds q* and g* 
depend only on the u-intersection (i.e. —) of the line cxi — cum — Ci 2 n = 0 and 
n-intersection (i.e. of the line <72 — C 21 n — C 22 n = 0 as well as the n-intersection 
(i.e. of the line ai — Cuu — C 12 V = 0 and u-intersection (i.e. of the line 
(72 — C 21 M — C 22 V = 0, respectively. We note that e 2 = (^)O) and 63 = (0, 
represent the two competitively exclusive states. When a = (3 = 1, the above 
observation together with the fact that the estimate of au{x) + /3v{x) given in 
Theorem |2 .1 1 does not explicitly depend on the propagating speed 9 gives a possible 
answer to Ql. 


(ii) Letting o; = 1 and /I —?• 0 in Theorem 2.1, we obtain 


0 < ■u(x) < max 


(7i (72 

LCii ’ C21 J 


max 


Letting (3 = 1 and o; —?• 0 in Theorem 2.1, we obtain 


0 < v{x) < max 


(72 (7i 

C 22 ’ C 12 


max 


di (^2 
-d2 di- 

d \ (^2 
d2 di 


( 2 . 6 ) 


(2.7) 
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(iii) When [S]: ^ ^ and ^ ^ holds, q* and g* are given by 



^2 0*^1 


' di d2 

min 

a —, p — 

min 

7^7 


C 21 Ci 2 _ 


-a2 di- 


and 


When fWl: ^ 

I- -1 cn C91 


max 

cri „ (72 
a —, p — 

max 

di 6^2 

757 


Cii C22_ 


.U2 di- 


and holds, q* and g* are given by 


and 



CTi 

p ^ 2 ’ 


'di 

(^2 

g* = min 

a —, 

B — 

min 


— 

Cii 

C 22 . 


-(^2 

di- 


0-2 

q ^ 1 ' 


'di 

(^2 

g = max 

cx , / 

C 21 

3 — 
C 12 

max 

-(^2 



(iv) By letting a = /? = 1, we have 

g* < u{x) + n(x) < g*, a; e M, 


where 




CTl 

(72' 


'(72 

(7l ' 


'di 

(^2 

g* = min 

min 

■Cll’ 

C21- 

, min 

-C 22 

C 12 - 

min 

■ d2 





' (7l (72' 


' (72 (7l' 


'di ^2 

max 

max 

1 

1 

0 1 

to 

1_ 

, max 

-C 22 ’ C 12 - 

max 

.(^2 ^1 - 


This answers Q2. 


(2.8) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


Let us return to the system of three species (1.1) and consider traveling wave solutions 

{u{y,t),v{y,t),w{y,t)) = {u{x),v{x),w{x)), x = y-et (2.15) 

satisfying the following boundary value problem 


di Uxx + 9Ux-\- u{ai — Ciiu — Ci2 n — Cis w) = 0, x G M, 

d2 Vxx + 9Vx + V ((72 - C 21 M - C 22 f - C 23 w) = 0, X G M, 


(2.16) 


4 Wxx+ 9 Wx+ W ( 0-3 - C 31 M - C 32 X - C 33 w) = 0, X G M, 
{u,v,w){-oo) = (^,0,0), {u,v,w){oo) = 


Here again 9 is the propagation speed of the traveling wave. As mentioned in the begin¬ 
ning of introduction, we will study the influence of an exotic species W on other native 
species U and V in terms of (2.16). The hrst question we shall ask is whether competitor- 
mediated coexistence occurs for u, v, and w in the system (2.16). If the three species do 
coexist under certain conditions, then what will be the prohles of u{x), v{x), and tc(x)? 
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The result in 


indicates that when w{x) is absent in (2.16), the system of two species 


(1.9) under certain conditions admits solutions {u{x),v{x)) having the prohles with u{x) 


being monotonically decreasing and v{x) being monotonically increasing. Moreover, we 
see from the prohles of u{x) and v{x) that u{x) and v{x) dominate the neighborhood of 
X = —oo and the neighborhood of x = cxd, respectively. These facts lead us to the expec¬ 
tation that, the prohle of w{x) must be pulse-like {w{x) is a pulse if w{—oc) = w{oc) = 0 
and w{x) > 0 for x G M) if it exists since w will prevail only when u and v are not 
dominant. 

To simplify the problem, we restrict ourself to the case of ai = cn in Section and 
denote a solution of ( 2.16| ) with cti = Cn by ( 1 , 0, 0 )-(m*, n*, 0) wave for convenience. 
Although we can hnd exact (1,0, 0)-(n*, n*, 0) waves for (2.16) (See Theorem 4.1 and we 
remark that when cti ^ Cn, a similar result remains valid.) under certain restrictions 
on the parameters, it remains an open problem, however, to establish the existence of 


solutions for (2.16) under more general conditions. In spite of this fact, when we consider 


the situation where the influence of the invading species W on the native species U and 
V is of no signihcance, i.e. Ci 3 ,C 23 ~ 0 in (2.16), the limiting case Ci 3 ,C 23 —)■ O’*" leads to 
the boundary value problem 


(2.17) 


di Uxx + 0ux+ u{ai — ciiu — Ci 2 v) = 0, x G M, 
d 2 Vxx + OVx + V {<^2 -C21U- C22 1 ^) = 0 , X G M, 

4 Wxx+ 0 Wx+w (0-3 - C31 M - C32 n - C33 w) = 0 , X G M, 

(M,n,'«;)(-oo) = (^,0,0), (u, n,'w;)(cx)) = (M*,n*,0). 

Under the assumption of the existence of solutions (n(x),n(x)) = ('u(x),h(x)) for the 


system of two species (1.9) 


), it will be proved in Section 5.1 that 


under certain conditions, a solution w{x) of the third equation in (2.17), i.e. the non- 
autonomous Fisher equation for w = w{x) 

dsWxx + Owx+ w{a 3 - C 31 U - C 32 V - C 33 W) = 0 , xgM (2-18) 

can be found applying the supers olution-sub solution method, thereby establishing the 


existence of solutions for (2.17). We remark that, as an application of the N-barrier 


maximum principle, upper and lower bounds of C 31 u -\- C 32 V are used in constructing 


supersolutions and subsolutions of (2.18). 


Theorem 2.2 (Existence of traveling wave solutions for three competing species) 

Assume either [S] or [W]. Suppose that there exist u = m(x) and v = h(x) which solve the 


first two equations in (2.17) for some 6 , di, ai, Ca {i = 1, 2), and Cij {i,j = 1, 2 with i 7 ^ j) 


and satisfy the boundary conditions {u,v){—oo) = (1,0), {u,v){oc) = {u*,v*). Let 




■ ( 7 i (72 ■ 


■ (72 ( 7 i ■ 


'di 6^2 

q = min 

C31 min 

■Cll ’ C21- 

, C32 min 

-C22 ’ C12- 

min 

.^2 di- 


and 

■ (Ti (72 
-Cll ’ C 21 

Assume that the following hypotheses hold: 


q = max 


C 31 max 


C 32 max 


(72 (7i 

C22 ’ C12 


max 


di d2 
d2 di 


(2.19) 


( 2 . 20 ) 
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[HI] 0-3 < C31, 0-3 < C31 u* + C32 V*; 

[H2] —C 33 K — q + as<0; 

[H3] 4 6*2 - 4 (C33 K + Qd 3 ) (-C33 K - 2(1^ - q + a^) < 0; 


[H41 K> K> 0 ; 


Then ( 2 . 17 ) has a positive solution {u{x),v{x),w{x)) with w{x) ^ 0 for a; G M and 
w{x) —)■ 0 as a; —)■ ±cx). Moreover, w{x) < w{x) < w{x) for x G M, where w{x) = 
K [1 — tanh^(x) ] and w{x) = K. 


Another consequence of the N-barrier maximum principle concerns the nonexistence 
of solutions for (2.16). In other words, we look for conditions on the parameters under 
which there exists no positive solution {u{x),v{x),w{x)) for (2.16). 


Theorem 2.3 (Nonexistence of traveling wave solutions for three competing 
species). Let Si = cxi C33 — (J3C13 and S 2 = cr 2 C33 — (73023. Assume that the following 
hypotheses hold: 


[Al] Si, S 2 > 0; 

[A2] either C21 Si > cn S2, C12 S2 > C22 Si or C21 Si < cn S2, C12 S2 < C22 Si; 


[A3] 


mm 


C 31 di min 


Si S 2 
Cll ’ C21 


, C 32 d 2 min 


S2 Si 
C22 ’ C12 


mm 


di (^2 
d2 ’ di 


> o's Css- 


Then (2.16) has no positive solution {u{x),v{x),w{x)). 


It will be clear from the proof in Section that [Al] and [A2] assure that the N- 
barrier maximum principle can be applied in proving Theorem |2.3[ We note in particular 
that when as is sufficiently small, [A3] in Theorem |2.3| clearly holds. From the viewpoint 
of ecology, the result of Theorem |2.3| states that as the birth rate as of the species W 
decreases below a threshold, the three species U, V and W no longer coexist. Intuitively, 
due to the weakness of the exotic species W, competitor-mediated coexistence cannot 
occur for the three species U, V and W in (2.16). 


3 N-barrier maximum principle: proof of Theorem |2.1| 


In this section, we use the notations p{x) = au{x) + (3v{x), q{x) = di au{x) -|- d 2 Pv{x), 
and F{u,v) = au {ai — cuu — cu v) + (3 v {a 2 — C 21 U — C 22 v) as in (1.13). We begin with 
a useful lemma. 


Lemma 3.1. For the quadratic curve F{u,v) = 0 in the uv-plane, 
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(i) ^ ^ o,nd ^ then F{u,v) = 0 represents a hyperbola; 

(ii) ^ ^ ^ then F{u,v) = 0 represents a hyperbola, a parabola, or 

an ellipse depending on the parameters in F{u,v) = 0. 


Proof. To prove (z), we calculate the discriminant V of the quadratic curve F{u,v) = 
au{l — u — aiv) + (3 kv {1 — a 2 U — v) = 0 to obtain 

F> = [a Ci 2 + 13 C21) — Act: 13 Cn 022- (3-1) 


Because of the assumption ^ ^ and ^ it follows that (a C12 +/d C2i)^ > 

4a/3C 12 C21 > AaI3 (C 22 cri (cn CT2) = Aa/3 Cn 022 - Thus > 0 and the quadratic 
curve F(u,v) = 0 is a hyperbola under the assumption and — > —. 

\ ^ y ^ cii C01 Coo C10 


C21 


C22 


C12 


For simplicity we let cxi = (72 = Cn = C22 = 1, C12 = and C21 = | to show (ii). 
Depending on the other parameters, it is shown in Fig 3.1 that F{u,v) = 0 represents a 
hyperbola, a parabola, or an ellipse in the uu-plane. 

□ 


Proof of Theorem 2^. An easy observation shows that, it suffices to prove for any a, (3 > 
0, we have 

< adi u{x) + (3 d 2 v(x) < g*, x eM., 


where 


and 


g* = mm 



■ ( 7 i (72 ■ 

, (3 d 2 min 

■ (72 ( 7 i ■ 


di (^2 

adi min 

■Cll ’ C21- 

-C22 ’ C12- 

min 

.d2 di- 
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■ (7i (72 ■ 

, [3 d 2 max 

■ (72 (7i ■ 


'di d2' 

max 

a di max 

■Cll C21. 

■C22 ’ C12- 

max 



(3.2) 

(3.3) 

(3.4) 


First of all, we prove (/) for the case of strong competition [SI: and — > —. 

_ ^ ir \ / Q SP >• J cii C21 C22 C12 

(3.5) 


Clearly, (3.3) in this case gives 
q{x) > min 


, 0-2 o j *^1 

a di —, [3 d2 — 
C21 C12 


mm 


d\ 6^2 


X G 


Vd2 diP 

The proof of the above inequality is divided into the following four cases: 


• \i d2> di, 


(z) when /3 Ui C21 ^2 > a <72012^1 and {Xi,X2, rj) 
Ai, X G M; 

{ii) when ( 3 aiC 2 id 2 < a <72 012^1 and (Ai,A2,z7) 
X G M. 


V C21 <12 ’ C21 ’ C21 02 / ’ ^ ^ ^ - 


/pjndi pui d2 
V C12 ’ C12 ’ 


/3o-i 

C12 


), q{x) > Ai, 


• It d2 < di, 

{Hi) when (3aiC2id2 > a <72 012^1 and (Ai, A 2 ,Z 7 ) = ^), q{x) > Ai, 

X G M; 
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iiv) 


when 13 ai C 21 d 2 < cta 2 C12 di and (Ai, A2, rj) = 
a: e M. 


\ ci2di ’ C12 ’ 


d2 

C12 di 


),q{x) > Ai, 


We first observe that the four cases can be reduced to the following two cases: 

• for (3 ai C21 d 2 >aa 2 C12 di, q{x) > min [^, , a; e M; 

• for P ai C21 d 2 <aa 2 C12 di, q{x) > min g], a; G M. 

Combining the two cases above leads to 


q{x) > min 


, 0-2 oj <^1 

a di —, [3 d2 — 
C21 C12 


mm 


di d2 


X G 


which is the desired result. The two inequalities in (2.4) give 


q"{x) + 0p\x) + F{u{x),v{x)) < 0, a: G 


(3.6) 


(3.7) 


For d 2 > di, we prove (i) by contradiction. Suppose that, contrary to our claim, there 
exists 2; G M such that q{z) < Ai and mina;g]Rg(a;) = q{z). Since -u,n G C^(M), we denote 
respectively by Z 2 and zi the hrst point intersecting the line adiu + (3 d 2 V = X 2 in the 
MU-plane, when the solution {u{x),v{x)) in the wn-plane travels from z towards 00 and 
—00 (as shown in Figure 3.2(a)). For the case where d < 0, we integrate (3.7) with 


respect to x from zi to z and obtain 


q'(z) — q'(zi) + 6 {p{z) — p{zi)) + / F(M(a;), n(a;)) da; < 0. 

J Zl 


(3.8) 


On the other hand we conclude: 


• since mina;g]R q(x) = q(z), it is easy to see q'(z) = adi u'(z) + /3 d 2 v'(z) = 0; 

• q(zi) = X 2 follows from the fact that zi is on the line adiu + f3 d 2 V = X 2 . On the 
other hand, when Zi moves a little towards 00 , q{zi+) is below the line adiu + 
(3 d 2 V = A2. This leads to q{zi + d) < A2 for any d > 0, and hence q\zi) < 0; 

• p{z) < p since 2; is below the line auX- (3 v = p] p{zi) > p since 2; is above the line 
a u + (3 V = p; 


it is readily seen that the quadratic curve F{u,v) = 0 passes throu gh t he points 
(0.0), (^,0).(0,g),and {u*,v*) in the wn-plane. Applying Lemma 3.1 1 it follows 
that F{u,v) = 0 is a hyperbola and F(u(x), v(x)) dx > 0 since F{u,v) < 0 as 
u,v > 0 are sufficiently large. 


Summarizing the above arguments, we obtain 

q'(z) - q'(zi) + 0 {p{z) - p{zi)) + [ F{u{x),v{x)) dx > 0, 


(3.9) 


fzi 


which contradicts (3.8). Therefore when d < 0, q{x) > Ai for a; G M. For the case where 
d > 0, integrating (3.7) with respect to x from ^ to Z 2 yields 

j‘Z2 

q {^ 2 ) - q {z) + 6 {j){z 2 ) - p{z)) + / F{u{x),v{x)) dx <Q. (3.10) 
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In a similar manner, it can be shown that q'{z 2 ) > 0, q'{z) = 0, p{z 2 ) > rj, p{z) < p, and 
F{u{x),v{x)) dx > 0. These together contradict (3.10). Consequently, (z) is proved 


for d 2 > di. For di = d 2 , we have q{x) = p{x) and (3.7) becomes 


p"{x) + 6p'{x) + F{u{x)^v{x)) xe 


(3.11) 


Moreover, when di = d 2 we take ^ ^ = h = i-®- the three lines adiu + 

(3 d 2 V = Ai, adiu + {3 d 2 V = A2, and au + (3v = p coincide. Analogously to the case 
of d 2 > di, we assume that there exists A G M such that p{z) < Ai and min2,gKp(x) = 
p{z). Due to mina,g]Rp(a;) = p(A), we have p'{z) = 0 and p'\z) > 0. By means of 
Lemma 3.1, F{u{z),v{z)) > 0. These together give p''{z) + 6p\z) + F{u{z),v{z)) > 0, 
which contradicts (3.11). Thus, p{x) > Ai for x G M when di = d 2 . As a result, the proof 
of (z) is completed. 

Clearly, we see from Figures 3,2(b), 3 j2(c) and 3,2(d) that the proofs of cases (zz), 
(zzz), and (zu) follow in a similar manner. This completes the proof of the case of strong 
competition [Si: ^ ^ and — > —. 

To give the proof for the case of weak competition [Wj: ^ <z ^ and — < —, which 

Q ir ^ ^ ^ Cii C21 C 22 Ci2 ’ 

leads to 

'di d2' 


q{x) > min 


, cn ni ^2 
a di —, (3 d2 — 
Cll C22 


mm 


-(^2 di 


X G 


(3.12) 


see 


by (3.3), we hrst see from Lemma [3T that F{u,v) = 0 in this case is a hyperbola 
for example Figure 3,1(a), Figure 3, l(b)| and Figure 3,1(f)) as in the case of strong 
competition [S], a parabola (see for example Figure 3 jl(c) and Figure 3,1(d)), or an 
ellipse (see for example Figure 3 1(e)) depending on the parameters in F{u,v) = 0. 


However, since we are only concerned with the curve F{u,v) = 0 in the hrst quadrant 

of the zzu-plane, it is readily seen from Figure 3.1| that for each of the three generic 

types of quadratic curves, we can construct an N-barrier for which the arguments used 

in proving the case of strong competition [S] remain valid. Moreover, in addition to the 

diffusion rates di, d2 and the coefficients a, /?, the lower bound of g(x) given by (3.5) 

under [S] is only involved with the minimal u-intercept of the lines ui — cuu — C12 v = 

0 and (72 — C21 zz — C 22 V = 0, i.e. u = —, and the minimal u-intercept of the lines 

cTi — Cll zi — C12 u = 0 and ct 2 — C21 zt — C22 z; = 0, i.e. v = —. Accordingly, for the case of 

^12 


weak competition [W], we conclude that (3.12) holds since the minimal zt-intercept of the 
lines (Ti — Cll zt — C12 z; = 0 and (T2 — C21 zz — C22 n = 0 is zz = and the minimal z;-intercept 
of the lines ui — cn zz — C12 z; = 0 and (72 — C21 zz — C22 n = 0 is z; = respectively. For the 
case of strong competition [S], (|3.5) is proved, whereas for the case of weak competition 
|W], we obtain (|3.12[). Combining the two inequalities (|3.5[) and (|3.12[) yields the lower 


bound of q{x) given by (3.3). This completes the proof of (J). 


As in the proof of (/), there are also four cases for the proof of {II) when the condition 
of strong competition [S] holds: 


If d2 > di, 

(z) when f3 (72 cn d 2 > ctai C 22 di and (Ai, X 2 ,p) = 
X G M; 

(zz) when /3 (72 cn d2 < aaiC 22 di and (Ai,A 2 ,z 7 ) = 
X G M. 


Pjlidi 2^) q(x) < Ai, 

V C22 «1 ’ C22 ’ C22 «! / ’ / — -*-5 
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• \i d2 < di, 


{Hi) when /3 ct 2 cii ^2 > aaiC 22 di and {Xi,X 2 ,r]) = ^), q{x) < Ai, 

X e M; 

{iv) when/3 as cii ds < acTiC22di and (Ai,A2,r7) = , q{x) < 

Ai, X G M. 


Combining the four cases above, it immediately follows that 

• for [3 (72 cii d 2 >aai C 22 di, q{x) < max for x G M; 

• for fd (72 Cii d 2 < aai C 22 di, q{x) < max for x G M, 

and hence we have 


g(x) < max 


1 *^2 

a di —, p d2 — 
Cll C 22 


max 


d\ d2 
-d2 d\ 


X G 


(3.13) 


Indeed, (3.13) holds, as is readily seen by employing similar arguments as above together 
with the N-barriers constructed in Figures 3,3(a), 3,3(b), 3 j3(c) and 3 |3(d) On the other 
hand, under the condition of weak competition [W], (3.4) leads to 


q{x) < max 


, (^2 o J *^1 

a di —, p d2 — 

C21 C12 


max 


d\ d2 
d2 d\ 


X G 


(3.14) 


which can be shown as in the proof of (3.12) by interchanging the roles of ^ (^, 
respectively) and ^ {■^, respectively) in (3.13). Therefore {II) of Theorem 


from (3.13) and (3.14). The proof of Theorem 2.1 is completed. 


2.1 


follows 


□ 


4 New exact (1, 0, i;*, 0) waves 


In this section, we always assume ai = Cn, unless otherwise stated. Looking for traveling 
wave solutions {u{x),v{x),w{x)) with the prohles of u{x) being decreasing in x, v{x) 
being increasing in x, and w{x) being a pulse (i.e., w{±oo) = 0 and w{x) > 0 for x G M) 
of (2.16) leads to the following ansatz ( 0 [3 lU |I9] ) for solving (2.16) 


'u{x) = i(M* + 1) + i(M* - 1) T(x), 

< n(x) = A:i(l + T(x))2, (4.1) 

^m;(x) = k 2 (1 -T2(x)), 


where T(x) = tanh(x), fci = ^ and k 2 is a positive constant to be determined. It 
is readily verihed that the ansatz ( 4.1[ ) satishes the boundary conditions at x = ±oo 
in (2.16). Since u{x), v{x) and w{x) in (4.1) are expressed in terms of polynomials in 
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tanh(a;) and ^ tanh(a;) = 1 — tanh^(a;), inserting (4.1) into the three equations in (2.16) 
gives 

diu^^ + eu^+u{ai-ciiu-ci 2 v-ci:iw) = Cio + CnT{x) + Ci 2 T^{x) + Ci3T^{x) , (4.2a) 


d 2 V:,:, + 9v^ + v{a2-C2iU-C22V-C23w) = v C 20 + C 21 T {x) + (22 T‘^ {x) , (4.2b) 


d 3 W^^ + 9w^ + w{a 3 -C 3 iU-C 32 V-C 33 W) = w Cso + C3i7'(a;) + C32 7'^(a:) • (4.2c) 

Equating the coefficients of powers of T{x) to zero yields a system of ten algebraic equa¬ 
tions: 


Cn = 0 , z = 0 ,l, 2 , 3 , 
C2* = 0, * = 0,1,2, 
C3* = 0, * = 0,1,2. 


It turns out that (4.2) can be solved to give 

Cll = O'!, 


Cl 2 = 


dial 


Cl 3 = 


ki {2di + 9) ’ 
2di9 + diai — Adf 


C 22 — 


C 23 — 


k 2 { 2 di 9 ) 

C 21 = 16^2 + 4:9 + a2, 

2di9 — 4(^2^ + diCT2 -|- 8did2 — 9'^ 
ki { 2 di + 9 ) ’ 

2di9 — 10 (^ 2 ^ + dia2 — 4did2 — 9"^ 
k 2 { 2 di 9 ) 

C31 = 4ci3 + 29 + CTg, 


C32 — 


diCTs -|- 4did3 — 9‘^ 


ki {2di + 9) ’ 

—6^30 -|- dia3 — Sdids — 9 '^ 
k2 {2di + 9) 

The result obtained is summarized in the following theorem. 


C33 — 


(4.3) 


(4.4a) 

(4.4b) 

(4.4c) 

(4.4d) 

(4.4e) 

(4.4f) 

(4.4g) 

(4.4h) 

(4.4i) 


14 



















Theorem 4.1 (Exact (1, 0, 0)-(m*, v*, 0) waves). The boundary value problem (2.16) 


with (Ti = cii admits an exact solution of the form (4.1) provided that (4.4) holds 


Theorem 4.1 asserts that under certain conditions imposed on the parameters, i.e. 


under (4.4), we can hnd exact (1, 0, 0 )-(m*, n*, 0) waves of (2.16) and the exact waves are 


polynomials in tanh(a:). To illustrate Theorem 4.1, let us choose ki = k 2 = di = d 2 = 


ds = 1, 9 = 3, and ai = a 2 = cr^ = 41 in (4.4). This gives Cn = 41,ci2 = ^,Ci3 = 


|1,C21 = 69 ,C 22 = f,C 23 = |,C31 = 51 , C32 ^^,033 = | and {u*,v*) = (|, 4 ). The 
resulting exact (1, 0, 0 )-(m*, n*, 0) wave is given by 


'u{x) = § ~ 5 tanh(a;), 
v{x) = (1 + tanh(a;))^, 
^w{x) = 1 — tanh^(x). 


(4.5) 


The prohles of u{x), n(x), and w{x) are shown in Figure 4.1 We conclude this section with 


in [19], where the ansatz (4.1) for solutions is 


the remark that the ansatz (4.1) for solutions of (2.16) is inspired by the one proposed 


u{x) = 2 (m* + 1 ) + - 1) T{x), 

v{x) = ki{l + T{x)Y, 


(4.6) 


when w is absent in (2.16) 


5 Applications of the N-barrier maximum principle 


5.1 Application to the existence of three species traveling waves: 


proof of Theorem 2.2 


To prove Theorem 2.2, we hrst observe that the third equation in (2.17) can be regarded 


as a non-autonomous Fisher equation when u = u{x) and v = v{x) are given, say u = u{x) 
and V = v{x). Then the non-autonomous Fisher equation (in, m, EH) is 


ds w^^ + Ow^ + w (0-3 - C31 M - C32 n - C33 w) = 0, x e 


(5.1) 


In order to hnd a solution of (5.1), an approach based on the supersolution-subsolution 


method is employed. To this end, we introduce supersolutions and subsolutions, w = 
w{x) is said to be a supersolution of ( |5.1 ) if it satishes the differential inequality 

ds Wxx + 9w^ + IV (0-3 - C31M - C32 h - C33 tc) < 0, X e M. (5.2) 


Similarly, a subsolution w = w{x) is dehned by reversing the inequality in (5.2). The 


following lemma is helpful in constructing non-trivial solutions of (5.1). 
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Lemma 5.1 (|5)- Suppose that w{x) is a hounded solution of 

dsWxx + Ow^ + w (p{w,x) = 0 , X G M, (5.3) 

where (^( 0 , x) — )■ —as x ^ ±oo for some constants > 0 . Then w{x) 0 as 

X —)■ ±oo. 


To construct a pair of subsolution and supersolution of (5.1), we employ the tanh 
method. 


Proof of Theorem \2.2 . Due to [H4], we clearly have w{x) < w{x) for x G M. First of all, 
we show that w{x) and w{x) are a subsolution and a supersolution of (5.1) respectively. 
Indeed, a straightforward computation gives us 


d2.Wxx + 0w^ + w{a^- C31 U-C32V- C33 w) 

= w [ — C33 K — 2 ds — q + as — 2 9 tanh(x) + (C33 K + 6ds) tanh^(x)] > 0, 

d3Wxx + dw^ + w{as- C31 U-CS2V- C33 w) 

= w [-cssK - q + as) <0 (5.4) 


under hypotheses [H2] and [H3]. By means of Theorem 2.1, we have used in the last 
two inequalities an estimate of C3i'u(x) + Cs 2 v{x), i.e. 


Q < C31 u{x) + C32 h(x) < g, X G M. 


(6.5) 


The existence of a solution w = w{x) for ( 5.1| ) lying between the subsolution w{x) and 
the supersolution w{x) constructed above follows from Theorem 2.8 in [2S]- In view of 


[HI], we hnally employ Lemma 5.1 to conclude that the solution w{x) of (5.1) has the 
asymptotic behavior hm|a;|^oo w{x) = 0. This completes the proof. □ 


5.2 Application to the nonexistence of three species traveling 


waves: proof of Theorem 2.3 


Proof of Theorem 2.3[ We prove Theorem 2.3| by contradiction. Suppose that, to the 
contrary, there exist u{x),v{x),w{x) > 0 , x G M satisfying (2.16). Since w{x) > 0 for 
X G M and w{±oo) = 0 , there exists xq G M such that max^jg® t(;(x) = w(xo) > 0 , 
w"(xo) < 0 , and w'{xo) = 0 . Due to ds Wxx + 6wx + w{as — csiu — C32 v — C33 w) = 0 , we 
obtain 

(^3 - C31 u{xq) - C32 n(xo) - C33 tc(xo) > 0 , ( 5 . 6 ) 


and hence 


w(x) < w(xo) < -(cr3 - C31 m(xo) - C32 n(xo)) < —, X G 

C33 C33 


(5.7) 
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As a result, we have 


di Uxx + Oux + u{ai — Ci3 (T3 — cii M — C12 u) < 0 , a; G M, 

d2 +0V:2 + v(/T2- C23 CT3 €33^ - C21 U - C22 v) < 0 , X £ R. 


(6.8) 


Because of [Al] and [A2], we can apply (/) of Theorem 2.1 to (5.8). Indeed, [Al] 
assures the positivity of Ui — and a 2 — 023 ( 730 ^^, whereas the assumption of 

strong competition [S] or the assumption of weak competition [W] for the nonlinearity 
in (5.8) follows from [A2]. Consequently, (/) of Theorem 2.1 gives ns a lower bound of 
C 3 iM(a:) + C 32 v{x), i.e. for x G M, 



'El S2' 

, C 32 d 2 min 

'^2 El" 


'di ^2' 

C 31 di min 

■Cii ’ C 21 - 

-C 22 ’ C 12 - 

min 

-d2 di- 


C 31 u{x) + C 32 v{x) > C33 min 

The condition [A3] then yields 

C31 u{x) + C 32 v(x) > 0-3, X G M, 


(5.9) 

(5.10) 


which contradicts (5.6). This completes the proof of the theorem. 


□ 
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Figure 3 . 1 : Red line: ai — cn u — C12 u = 0; blue line: (T2 — C21 u — C22 v = 0; green curve: 


au(ai-cnu-ci 2 v)-t-/ 3 v (o^C2iu-C22v) = 0. iTi = cr2 = cii = C22 = l,ci2 = 5,021 = (a) 


a = 5, /3 = 4 (hyperbola). 

(b) 

a = 2, /3 = ^ (hyperbola). 

(c) 

a = 

(parabola). 

(d) 

a = 2, /3 

= ^ — 3\/6 ~ 0.1515 (parabola). 

(e) 

zooming out of (a). 
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(a) 


(b) 




(c) 


(d) 


Figure 3.2: Red line: ai — cn u — ci2U = 0; blue line: (T2 — C21 u — C22V = 0; green curve: 
a u (cJi — cii M — C12 v) + ( 3 v {a2 — C21 n — C22 v) = 0 ; magenta line (above): adiu + (3 d2V = X2; 
magenta line (below): adiu + f 3d2V = Ai; yellow line: au+( 3 v = r/; dashed curve: (u(x),u(x)). 
di = (ji = (T2 = cii = C22 = 1- (a) C12 = 2, C21 = 3, a = 17, /? = 18, and ^2 = 2 give Ai = 

A 2 = and T] = (b) C 12 = 2, C 21 = 3, a = 17, = 5, and ^2 = 2 give Ai = |, A 2 = 5, and 


r/ = 


17 


(c) C12 = 2 , C21 = 3 , a = 17 , = 18 , and ^2 = | give Ai = ^, A2 = -^, and V = 3 


C12 = 2 , C21 = 3 , a = 17 , /3 = 18 , and ^2 = 5 give Ai = |, A2 = |, and d = §• 


(d) 
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(a) 


(b) 




(c) 


(d) 


Figure 3.3: Red line: cti — chu — C 12 n = 0; blue line: (T 2 — C 21 n — C 22 u = 0; green curve: 
a u (cJi — cii n — C 12 v) + (3v {a 2 — C 21 n — C 22 v) = 0 ; magenta line (below): adiu + (3d2V = A 2 ; 
magenta line (above): adiu + f3d2V = Ai; yellow line: au+j3v = r/; dashed curve: {u{x),v{x)). 


di = cTi = (T 2 = cii =£22 = 1- [^ci 2 = 2, C 21 = 3, a = 17, /3 = 18, and ^2 = 2 give Ai = 72, 
A 2 = 36, and ry = 36. |(b)| ci 2 = 2, C 21 = 3, a = 17, /? = 5, and ^2 = 2 give Ai = 34, A 2 = 17, and 
ry = 17. (c) C 12 = 2, C 21 = 3, a = 17, /3 = 33, and ^2 = | give Ai = 33, A 2 = 22, and ry = 33. 
(d) C 12 = 2 , C 21 = 3, a = 17, /3 = 18, and ^2 = 5 give Ai = 34, A 2 = 17, and ry = 34. 
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Figure 4.1: Profiles of the solution {u{x),v{x),w{x)). 
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